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'nI" ; Abstract 
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^ , We present an analytical approximation for nonlinear dynamics of trapped Bose- 

O 
D ' condensed gases. The new approximation is a substantial improvement over the 



Thomas-Fermi approximation and is shown to be applicable for systems with 
a rather small number of atoms A^. The calculated aspect ratios after ballis- 
tic expansion are found to be in good agreement with those observed in recent 



Q I experiments. 
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The newly created Bose-Einstein condensates (BEC) of weakly interacting alkali-metal 
atoms [1] stimulated a large number of theoretical investigations (see recent reviews [2]). A 
mean-field approach, based on the time-dependent Gross-Pitaevskii (GP) equation [3], is the 
most widely used theory for nonlinear dynamics of trapped Bose-condensed gases at zero 
temperature. 

In the limit of a large number of atoms N, the determination of the condensate wave 
function is simplified by neglecting the kinetic energy term. This approximation is known 
as the Thomas- Fermi (TF) approximation [4,5]. It has been used quite extensively [7-13], 
including the explicit time evolution of the condensate ( shape of profiles, aspect ratio etc.) 
during the expansion after switching off the trap [7,9]. We note that the validity of the TF 
approximation depends not only on the number of atoms A^, but also depends on properties 
of the traps. 

In this letter we develop an analytical approximation which is a substantial improvement 
over the TF results for rather small number of atoms A^. The aspect ratios calculated from 
the new approximation method are also found to be in good agreement with those observed 
in recent experiments [14]. 

In the mean-field approximation, the ground state energy of the system is given by the 
Ginzburg-Pitaevskii-Gross (GPG) energy functional [3,15] 
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<^\Y^H,\^>+^J\^\Uf, (1) 



«=1 
with 

and g = Airfi^a/m, where a is the S'-wave scattering length, A^ is the number of atoms in the 
BEC, and \E' is the condensate wave function normalized as / | \E' p (ir*= 1. 
We introduce an auxiliary Hamiltonians 

^. = ^v^+^(l-7.K^ (3) 

where 7j are parameters, < 7j < 1, and rewrite Eq.(l) as 

I =< M/ I J2{H,-H,) \^> + <^\J2H,\^>+^J\^ r dr. (4) 
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Omission of < \1/ | YA=i{Hi — Hi) | ^ > in Eq.(4) yields our approximation for the ground 

state, 

f^ = ^^Vli + fJ'TF{i^x,^y,U;z), (6) 

i=l 

and 

p{r) = p'^^{u^,Uy,u;,,r), (7) 

where cDj = ^/l — 'jiUJi. Etf^I^tf^ and p^^ are the Thomas-Fermi energy, chemical potential, 
and density, respectively, which are given by 
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ETF{oOx,^y,uJ^) = —N[—uj^u;yUJ^m^/^gN]'^/\ (8) 
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HTF{^x,^y,^z) = -[—UJx^y^zm^^'^ gNf'^ , (9) 

and 

p^^{u.^uy,u.,r) = ^^^^^!;^"^-V -E( J^n^a-El^F)^), (10) 

^9 i=i ^i i=i ^i 



where 



IUTF\2 ^ '2f^TF{^x,^y,^z) ,^^. 

* mujf 



Projecting | ^ > on the complete basis states \ n >, obtained from hi = ~ ^ f? + ™'^' ^' x'i ? 

i 

and hi \ n >= e„ | n >, we get 

< ^ I /i . I ^ >= E^« |< ^ I n >< n I ^ >| > ei = ^'^ . (12) 

n ^ 

Therefore, we conclude that our approximation for energy, given by Eq.(5), is a lower bound 
to the ground state energy, Eq.(4). Therefore a set of the optimal values of parameters 7j 
which maximizes the energy, Eq.(5), will yield an optimal value for the ground-state energy 
given by 

J = l 4 = 1 

Since the TF approximation corresponds to the case of 7j = 0, we have 

EtF < E < Eexact: (14) 

where Eexact is the exact mean-field energy. 



To study the validity of our approximation, we consider an example of the ground state 
of ^"^Rb atoms in a harmonic trap, as investigated in Ref.[16] with the S-wave triplet-spin 
scattering length a = lOOas, where as is the Bohr radius, the axial frequency ujz/27i = 
220Hz, and assymetry parameter A = oj^I^l = V8, where to^ = ujy = ^±- 

Using our approximation, we calculate the energy per particle, E/N, the chemical poten- 
tial /i, and the average transverse sizes vx'^ and vertical sizesv z^, using Eqs. (5-6,13). The 
calculated results are compared with those obtained from the numerical solutions of the GP 



equation, Enum/N, fXnum, yx^nurn and y z^num [16] in Table I, and with those obtained in 



the TF approximation Etf/N, htf, yx^TF and y z^tf in Table II. These comparisons show 
that our analytical approximation greatly improves the TF results for a rather small number 
N . For 100 < N < 20000, the difference between our results and those of the numerical 
solution of the GP equation [16] are less than 3%. 

Let us now turn to application of our approximation for the time-dependent problems. 
Consider the BEG with the time-dependent harmonic potential Vt = {rn/2)Y^^=iOjf{t)xf. 

In Ref.[7] the following anzatz 

^(r, t) = exp[-t(3{t) + «m 2^ — -— ] (15) 

V A1A2A3 j^^i zn Ai 

has been used for the solution of the time- dependent GP equation 

ihj- = A^ + \/i(f,t)^ + iV^ I ^ |2^ (16) 

(a similar treatment for isotropic traps has been developed in Ref.[17]). We choose (3 and A 
satisfying the following equations 

^/3 = E ^^ - ^^^f??'^^\ m) = 0, (17) 

fc=i 2A;, A1A2A3 



and 



A.A. + Xlu^m - ^ = ^\ ^f\ A.(0) = 1, A.(0) = 0, 
Afc A1A2A3 



with Uk = uJk{0), and uJk = y^T^^^yj^uOk- For a special case of 7^ = 0, the above choices of (3 
and Afc reduce to those of Ref.[7]. With substitutions of Eqs. (15,17-18), Eq.(16) becomes 



^^d^ = Y. ^^V^$ + ^^[-/iTF(cJ.,cJ„^.) + ^ E^'^fc + Ng I $ 11$, (19) 



with the initial condition $(r*, 0) = \l/(r*, 0), where \l/(r, 0) is a solution of the time- independent 
mean-field equation 

-— A^(f, 0) + Vt{f, 0)^(f, 0) + A^^ I ^(f, 0) |2 ^(f, 0) = ii^{f, 0). (20) 

By neglecting Yj\=i ^''7'^^'° in Eq. (19), we obtain a generalization of our approximation, 

k 

Eqs.(5-7,13) to the time-dependent problem 

I ^J.{^f\ |2_ PTF{uJx,i^y,^z,Xi/Xi{t),X2/X2{t),X3/X3{t)) 

' ^''^'" Ai(t)A2(t)A3(t) ' ^''' 

where all the dynamics is in the evolution of the scaling parameters A^, Eq.(18). 

For the case u^ = ujy = ujj_ and Ai = A2 = A_l, the aspect ratio of the cloud in our 
approximation is given by 
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As an example, we consider application of the above results to the experimental data 

with "^^Na atoms obtained in the loffe-Pritchard type magnetic trap with radial and axial 

trapping frequences of c<jj_/(27r) = 360 Hz and ujz/{2n) = 3.5 Hz [14], respectively. In our 

analysis we use a = 2.75nm, t = 4 ms, and a/a± = 2.488 x 10"^, where a± = Jh/muj±. As 

in Ref. [7] we consider a sudden and total opening of the trap at t = 0. For this case Eq. 

(17) becomes 

d^\± 7_L , 1 - 7± 



dr^ Ai AiA, 



(23) 



where r = u±t and e = uJz/uj± ^ 1. 

To the zero-th order in e, we have A^ = 1 and A_l = a/1 + t^- For the experimental 
conditions [14], the terms in e^ are negligible. 

In Table HI, we give the calculated values of the aspect ratio R{t) of the "^^Na atoms 
cloud, after ballistic expansion of t = 4 ms, and the calculated values of parameters 7_l and 7^ 
which we fix from Eq.(13) , with Ui = uji{0). One can easily see that the TF approximation 
is valid (7^ ~ 0) along the long axis of the cloud, but not in the radial direction (7^ 7^ 0), 
as has been noted already in Ref. [18]. 



Our calculated results for R{t) are compared with the recent experimental data [14] in 
Figure 1. This comparison shows that our predictions for the aspect ratio R{t) are in good 
agreement with experimental data obtained by the MIT group [14]. 

In conclusion, we have developed an analytical approximation which provides a sub- 
stantial improvement over the TF approximation for nonlinear dynamics of trapped Bose- 
condensed gases for a rather small number of atoms A^. The approximation is very useful 
since it provides an easy quantitive tool for the analysis of experiments on trapped condensed 
gases. 

We are grateful to the group at MIT for providing us with the experimental data. 



Table I. Comparison of the results of our approximation for the ground state of ^^Rb 
atoms, calculated from Eqs. (5-6,13) and those obtained from the numerical solutions of the 
GP equation [16]. Chemical potential and energy are in units of hu±, and length is in units 



of Jh/mLJ± 





E/N 


/i 


Va;2 


V22 


-C'riMm/-' ' 


f^num 






N 




\/z^ 

V num 


100 


2.63 


2.82 


0.78 


0.43 


2.66 


2.88 


0.79 


0.44 


200 


2.80 


3.13 


0.83 


0.45 


2.86 


3.21 


0.85 


0.45 


500 


3.22 


3.82 


0.94 


0.47 


3.30 


3.94 


0.96 


0.47 


2000 


4.49 


5.49 


1.21 


0.53 


4.61 


5.93 


1.23 


0.53 


5000 


5.99 


8.00 


1.46 


0.59 


6.12 


8.14 


1.47 


0.59 


10000 


7.63 


10.4 


1.68 


0.65 


7.76 


10.5 


1.69 


0.65 


15000 


8.84 


12.1 


1.83 


0.69 


8.98 


12.2 


1.84 


0.70 


20000 


9.84 


13.5 


1.94 


0.72 


9.98 


13.7 


1.94 


0.73 



Table II. Results of the TF approximation for the same case as Table I. 



N 


Etf/N 


fiTF 


\/x^TF 


^4f 


100 


1.44 


1.60 


0.68 


0.24 


200 


1.51 


2.11 


0.78 


0.27 


500 


2.18 


3.05 


0.93 


0.33 


2000 


3.79 


5.31 


1.23 


0.44 


5000 


5.47 


7.66 


1.48 


0.52 


10000 


7.22 


10.1 


1.70 


0.60 


15000 


8.49 


11.9 


1.84 


0.65 


20000 


9.53 


13.3 


1.95 


0.69 



Table III. Calculated aspect ratio of the '^^Na atoms cloud, using Eq.(22), after a ballistic 
expansion of t = 4 ms, as a function of A^. A^ is in units of 10^. 



N 


R{t) 


7± 


Iz 


1.2 


0.110 


0.354 


0.0 


0.8 


0.117 


0.426 


0.0 


0.4 


0.132 


0.551 


0.0 


0.1 


0.183 


0.766 


0.0 
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FIG. 1. Aspect ratio R of the "^^Na atom cloud after a ballistic expansion of t = 4ms, 
as a function of the number of atoms A^, with C(j_|_(0) = 27r x 360Hz, 072(0) = 27rx 3.5 Hz. 
Diamonds, dashed line, and circular dots represent the results of theoretical calculations using 
Eq.(22), the TF approximation, and experimental data from MIT group [14], respectively. 
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